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^ ' Abstract 
^ ■ 

Q^' The generalization of the hard thermal loop effective theory to anisotropic plasmas is described 

^. with a detailed discussion of anisotropic dispersion laws and plasma instabilities. The numerical 

QJ I results obtained in real-time lattice simulations of the hard loop effective theory are reviewed, 



both for the stationary anisotropic case and for a quark-gluon plasma undergoing boost-invariant 
expansion. 
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}£2 1 Introduction 
o 



1.1 Weakly vs. strongly coupled quark-gluon plasma 

00 ' The wealth of data harvested at the Relativistic Heavy Ion Collider (RHIC) [T] has led to a shift of 
paradigm in thinking about the quark-gluon plasma. The strong collectivity that is being observed, in 
particular in elliptic flow and jet quenching, is widely taken as pointing to a strongly coupled plasma 



^ which is qualitatively and quantitatively different from a parton plasma that can be described by 
perturbative quantum chromodynamics (QCD); for an opposing point of view see however Ref. [2]. 

In fact, already before the mounting experimental evidence for a strongly coupled quark gluon plasma 
(sQGP) at RHIC, perturbative QCD at finite temperature was in some difficulty describing lattice data 
on the thermodynamics of deconfined QCD. Lattice results for the thermodynamic pressure of QCD 
typically lead to a rather sudden rise of pressure and entropy to about 15-20 % of the Stefan-Boltzmann 
result for an ultrarelativistic gas of quarks and gluons, at a few times the transition temperature 
Tc. A straightforward first-order calculation in as in fact gives just this ballpark of deviations from 
the interaction- free result. Higher-order calculations require resummation of collective effects such as 
Debye screening and have been carried through to order a^lna^ O HI [5], but at face value they show 
hopelessly poor convergence of perturbation theory at all temperatures of practical interest (in fact up 
to ridiculously high temperatures ~ lO^Tc). 

However, it is now understood that the poor convergence properties of thermal perturbation theory 
beyond first order perturbation theory is at least to a large part signalling the need for more complete 
resummations of screening phenomena, since similar problems appear already in the rather trivial case 
of scalar 0(iV — > oo) models, which can be solved exactly and where the only effect is the generation of a 
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Figure 1: Comparison of the perturbative 3-loop result for the pressure of pure-glue QCD 
with lattice results from [6] (thick dark line; thickness representing roughly the statistical 
error). The light-gray band represents the result to order when the renormalization scale 
is varied between ttT and 47rT and the medium-gray band shows the reduced renormaliza- 
tion scale ambiguity when higher order terms are resummed through the parameters of the 
effective field theory of dimensional reduction. The lines marked "PMS" (for principle of 
minimal sensitivity) and "FAC" (fastest apparent convergence) correspond to two optimiza- 
tion prescriptions for the renormalization scale. [7j 



thermal mass [S] . Indeed, already a minimal resummation of the Debye mass beyond strict perturbation 
theory together with a simple optimizations of the (huge) renormalization scale dependence \l\ (see 
Fig. d]) gives a remarkably good description of the (continuum-extrapolated) lattice results down to 
about 2.5Tc. Fig. [2^ shows that a description of the entropy of (pure glue) QCD in terms of hard-thermal- 
loop (HTL) resummed quasiparticle propagators together with a standard 2-loop running coupling 
(with renormalization scale varied about the Matsubara scale 27rT) also gives a good description of the 
thermodynamics above about ST^ [9j (see also Ref. [lOj and references therein) |1| 

Meanwhile, new techniques have become available that allow the analytical treatment of strongly 
coupled gauge theories. In particular. A/" = 4 supersymmetric Yang-Mills theories (SYM) in the limit 
of large number of colors and strong 't Hooft coupling A = g^N = AnasN is now often taken as a 
model for hot QCD. On the basis of the AdS/CFT conjecture [161 113 one can make predictions for the 
otherwise inaccessible strong-coupling regime, notably for real-time quantities such as the specific shear 
viscosity r]/S. While ri/S ~ (A^logA)"^ 1 at weak coupling [HI [19], the AdS/CFT correspondence 
gives r]/S = l/An + 0(A"'^/^) at large 't Hooft coupling [2ni[2I], and such (extremely) low values seem 
to be indeed required in viscous hydrodynamics models of heavy-ion collisions at RHIC [22] 

For A/" = 4 SYM at finite temperature, no lattice results for the thermodynamic potential are 
available that would allow one to test either weak coupling results or the (conjectured) strong coupling 
results at finite coupling (which does not run because of conformal invariance). However, successive 
Fade approximants that interpolate between the known weak and strong coupling results give smooth 
and seemingly robust extrapolations, as shown in Fig. [2]d [23]. A comparison of the QCD results for the 
entropy with those for A/" = 4 SYM suggests that the strong-coupling expansion is no longer working 
well when the deviations from the Stefan-Boltzmann result (So) are less than some 15%, as is the case 
for QCD at temperatures above ST^. Such temperatures are expected to be reached in the heavy- 



^HTL quasiparticle models have also been used as phenomenological models down to the phase transition by adding 
fitting parameters in the running coupling [111 1121 113) , sometimes extended by incorporating quasiparticle damping in a 
form motivated by HTL perturbation theory [TH [TS] . 
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Figure 2: The left panel (a) compares the lattice result for the entropy of pure-glue QCD 
of [B] with a next-to-leading HTL approximation (NLA) where the quasi-particle entropy 
of HTL quasiparticles is corrected by next-to-leading order corrections for the asymptotic 
thermal masses in a simple model for a gap equation (the parameter c\ is introduced to 
restrict this correction to hard modes above the scale a/27iTV?2dCa) • The upper and lower 
dashed lines correspond to renormalization scales flus = 47rT and ttT, respectively. The 
right panel (b) shows the weak and strong coupling results for the entropy density of A/" = 4 
SYM theory together with the NLA results obtained in analogy to QCD, but as a function 
of A, which here does not run. The dashed and full heavy gray lines represent the Fade 
approximants and -R[4,4] which interpolate between weak and strong coupling results to 
leading and next-to-leading orders, respectively [23] . 



ion experiments at the LHC, which may eventually get to probe specific perturbative features of the 
quark-gluon plasma. But already in the case of RHIC physics, it is clearly mandatory to improve our 
understanding of both the weak and the strong coupling asymptotics of the quark-gluon plasma, as the 
truth will most probably be somewhere in between, perhaps quite far from either. 

In some respect, a weakly coupled quark-gluon plasma (wQGF) can be more difficult to describe than 
a strongly coupled one, in which strong coupling dynamics wipes out any structures from quasiparticle 
dynamics. In a wQGF, a small coupling leads to a hierarchy of scales, which need to be disentangled by 
assuming <^ 1. Of course, for quantitative predictions, extrapolations to ~ 1 have to be considered, 
but this is no more or less problematic than the equally bold extrapolations theoretician need to do by 
starting out from infinite ('t Hooft) coupling. 



1.2 Scales of wQGP 

In thermal equilibrium, the primary scale is the temperature T, which determines the mean energy of 
( "hard" ) particles in the plasma. The "soft" scale gT with -C 1 is the scale of thermal masses that 
determine the plasma frequency, below which no propagating modes exist, the Debye screening mass, 
and the scale of Landau damping. In Feynman diagrams, these effects come from one-loop contributions 
with the highest power of the loop momentum, cut off by the thermal distribution function, and are 
therefore called hard thermal loops (HTL). To determine the physics at softer scales, one typically needs 
to consider the effective HTL theory and its resummed propagators and vertices. 

At the scale g^T one encounters in fact a barrier for perturbation theory when the magnetostatic 
sector is involved, since the latter is characterized by a completely nonperturbative dimensionally re- 
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duced Yang-Mills theory, which by itself can be (comparatively easily) handled by lattice methods (the 
difficult part being the perturbative matching to the four-dimensional theory). Quantities such as the 
color relaxation or gluon damping rate f2M ES] and also the Debye screening length at next-to-leading 
order [26] involves logarithms of the coupling, whose coefficient can be calculated perturbatively, but 
which are nonperturbative beyond the leading log. 

The scale g^T is the one characteristic of large-angle scattering and also of the inverse shear viscosity 
T'^/rj, whose calculation at weak coupling [HI [19] requires resummations beyond, but building upon, 
hard thermal loops. 

When it comes to nonequilibrium physics, it turns out that the parametrically most important 
phenomena are plasma instabilities [271 [28] that appear already in the coUisionless limit, on the level of 
the HTL masses (and only those of gauge bosons [29]). As pointed out by Arnold et al. [30], this leads 
for instance to the necessity of a complete revision (still to be worked out) for systematic perturbative 
scenarios of thermalization ^SjJ . Intriguingly, plasma instabilities may also be responsible for anomalous 
contributions [22] to the inverse shear viscosity whose standard weak coupling result appears much too 
small (which is one of the key arguments in favor of sQGP). 



2 Hard thermal and hard anisotropic loops 

2.1 Hard (thermal) loop effective theory 

The effective HTL theory of QCD [531 El] is given by the collection of all one-loop diagrams which are 
proportional to T^, assuming that all external momenta are soft and therefore negligible compared to 
T. For external momenta ~ gT, the HTL self energy and vertex diagrams are parametrically of the 
same order as the corresponding tree-level quantities and therefore need to be resummed completely. 
A remarkably compact form of the HTL effective action can be written down formally [35j, which for 
gluons reads 



>^^HTL(a;) = 




Here ~ gT is the Debye mass in the HTL approximation, (■ ■ ■)v a normalized average over the 
directions vinv'^ = (1, v) with = 1, and Z)^ the covariant derivative in the adjoint representation. In 
the non-Abelian case, this leads to an effective action which is nonlocal and nonpolynomial, containing 
an infinity of vertex functions. The reason for nonlocality is that, in contrast to other examples of an 
effective field theory, one is integrating out stable real particles rather than virtual ones. Albeit elegant 
in form, the effective action ([T]) is not well-defined as it stands, because it still requires boundary 
conditions to be taken into account after the extraction of self energy and vertex diagrams. 
Alternatively, the HTL effective theory can be derived from the effective field equation 

D,Fr = m = g J ^,^,SUp, X, t) , (2) 

where 5 fa is a linearized perturbation of the color-neutral background (thermal) distribution of color- 
carrying hard particles v ■ 9/o(p,x, t) = 0, = p^/p^, satisfying gauge covariant Boltzmann-Vlasov 
equations [36j 

V ■ D 5/„(p, X, t) = gv^Frdlf^foip, x, t) = -g^E, + v x B„) ■ Vp/o . (3) 

HTL vertex functions are then defined as the functional derivatives of the induced current j[A] with 
respect to A^. In particular, the HTL gauge boson self energy is given by the linear terms in yl°, which 
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Figure 3: Location of poles in A^^ and A^,. The right part with > corresponds to 
propagating normal modes, the left part to (dynamical) screening. 



in Fourier space and imposing retarded boundary conditions reads 

nk)U = 9' J ^,v%^fo{p) (^g,, - ^^^^ A\k) = W'{k)AAk) . (4) 

Written as gauge-covariant Boltzmann-Vlasov equations, the HTL effective theory allows us to 
immediately generalize to nonthermal cases. When the background distribution function /o is stationary 
and homogeneous, i.e. only dependent on momentum, one can in fact still write down an effective 
action similar to ([T]) [37j. The main difference to the thermal situation is that Vp/o need no longer 
be proportional to p and thus v. Whereas in the thermal (or just isotropic) case, the magnetic term 
in the Vlasov equation ([3]) drops out, in the anisotropic case magnetic interactions lead to much more 
complicated and rich dynamics. 

2.2 Isotropic dispersion laws 

In the isotropic case, the polarization tensor U^'^ in (j4]), which is transverse with respect to the 4- 
momentum k, contains two structure functions corresponding to the two symmetric and transverse 
tensors A^^, = g^^ - - B^^, B^^ = with tt^ = {g^^ - !^)u'' and plasma rest-frame velocity 

Ua = Ut = ^ vn^" = 1 (n% - n^) , = n,. = -^Hoo (5) 
n% = ml, noo = "^^1 - ^ In ■ (6) 

As a consequence, the gauge boson propagator, which in Landau gauge reads 

-G,, = AtA^, + AlB^,, At = [P-Ut]-\ AL = [k'-UL]-\ (7) 

has two branches of poles, corresponding to two different dispersion laws, A^^ = for spatially trans- 
verse polarizations, and A^^ = for polarizations along u. 

The two branches of poles are shown in Fig. [3], plotted in a;^ over k^ in order to show poles cor- 
responding to propagating modes (k^ > 0) and poles corresponding to screening (k^ < 0), which are 
analytically connected, on the same plot. Propagating modes are seen to exist for uj^ > u'^^ = ml/3, 
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Figure 4: The magnetostatic screening mass squared as a function of the anisotropy param- 
eter ^ as given in ([T^ . 



with the transverse mode approaching the mass hyperboloid /c^ = = m'j-,/2 asymptotically. The 
poles of branch L approach the light-cone exponentially, and in doing so, their residues vanish also 
exponentially, signifying a purely collective mode. For frequencies below the plasma frequency, u < cDpi, 
there is screening of both electric and magnetic fields, with screening lengths depending on frequency. 
In the static limit, mode L has inverse screening length m^), the so-called Debye mass. This corresponds 
to the screening of (chromo-)electric charges in the medium. By contrast, mode T has infinite screen- 
ing length in the limit of vanishing frequency, corresponding to the absence of (chromo-)magnetostatic 
screening in the HTL approximation. 

2.3 Anisotropic dispersion laws and plasma instabilities 

In the anisotropic case, the fact that 11'^'' is symmetric and II'^'^ fixed by transversality k^U^'^ = would 
lead to 6 structure functions in general. Assuming that there is just one direction of momentum space 
anisotropy, n = (0,0, 1) (i.e., axisymmetry around the 2;-axis), one can define 4 symmetric tensors for 
n*-', corresponding to 4 independent structure functions. 
Defining spatial tensors 

j^ij = _ k^y/e, = k'y/k^, = n'h^/h\ d'^ = k'h^ + yfi\ (s) 

with n* = A^^n^ and decomposing the spatial part of II^'^ according to 

n^^' = aA'^ + pB'^ + -iC'^ + 5D'^ (9) 

one finds that the gluon propagator in temporal axial gauge reads |38] 

A(fc) = AtA + (k^ - + a + 7)A£B + [{(3 - u'^)Ac - Ar]C - ^A^D (10) 
At(A;) = [k2-cu' + a]-\ ^c{l^) = [{\^ - uj^ + a + -i){(3 ~ uj^) -\<i^h^6^]-\ (11) 

This propagator contains one branch of poles from A^^^ = 0, and in general two branches from A^^ = 0, 
except when k || n and thus n = 0, in which case there are only two branches in total. 
A special important case for an axisymmetric distribution function is given by 

/(p) = /iso (p' + e(p ■ n)^) (12) 
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Figure 5: Illustration of the mechanism of filamentation (Weibel) instabilities (adapted 
from Ref. [39]). Charged plasma particles moving transversely to the wave vector of a 
seed magnetic field tend to separate in counterstreaming currents which in turn amplify the 
magnetic field. 



with some isotropic (for instance thermal) function /iso- The anisotropy parameter ^ is allowed to take 
values from —1 to oo, with — 1 < ^ < corresponding to prolate (cigar-shaped) momentum anisotropy, 
and < ^ < oo corresponding to oblate (squashed) distributions. The polarization tensor in the 
(anisotropic) hard loop (HL) approximation (j4]) can then be evaluated in closed form |38j|. Changing 
variables + ^(p ■ n) 



p^ gives 



m 



m 



dn j f ' + ^(v.n)r2' 



An 



-V 



9 



27r 



l + e(v.n)2)2 ^ 

app — — 

dp 



6'' + 



v^k'- 



k ■ V + ie 



(13) 



As a simple special case, let us consider structure function a{k) for the case that k || n and in the 



static limit where a{uj = 0) 
squared. One easily finds 



2^^ 



■0) 



mz 



has the interpretation of magnetostatic screening mass 



kiln 



-0(-e)-'/'atanh(-0^/^ 
Or^/^arctan^i/2 _ i] 



for ^ < 
for ^ > 



(14) 



which is plotted in Fig. |H 



For ^ 



0, 



vanishes, reproducing the result that there is no HTL magnetostatic screening 



mass. However, for ^ < (prolate momentum anisotropy) we obtain a nonvanishing HL magnetostatic 
screening mass, and for ^ > this mass turns out to be imaginary, signalling an instability. 

This magnetic instability can be identified with the so-called Weibel instability 00] and is illustrated 
in Fig. [5l In an extremely oblate momentum distribution, where all hard particles move in planes 
orthogonal to the small seed magnetic field with wave vector in the z-direction tends to sort 

streams of charged particles in a way which leads to an induced current that reinforces the initial 
field. The induced current and the magnetic field grow exponentially, until the magnetic field is large 
enough to bend the trajectories of the hard particles significantly, which in an Abelian plasma provides 
a mechanism for fast isotropization. Of course, in the HL approximation, we can only study the onset 
of such instabilities, because the backreaction on the background distribution /o is neglected. 
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Figure 6: Anisotropic dispersion laws for k || n and three values of ^: (isotropic), —0.5 
(prolate), and 5 (oblate). Full and dashed lines represent poles of and A£, respectively. 
The magnetic Weibel instability corresponds to the poles for negative uP' and positive k^ 
that appear for positive ^. 



In fact, there are also instabilities when ^ < 0, but to see those we have to consider wave vectors 
pointing away from the z-axis. For small ^, one finds that in the static limit there are the following 
poles for real wave vector k with /c^/lkl = cos 6', 



A^i 



: 
: 



kVm^ = -i cos2 + O(e') 
3 

kVm^ = -i cos2^ + 0(e2) 
3 



(15) 
(16) 



Thus, when ^ < 0, there are space-like poles in A^ when 7r/4 < Q < Stt/A, corresponding to electric 
(Buneman) instabilities [50] . 

Returning to the case of oblate distribution functions, > 0, and restricting again to k || n (i.e., 
I cos^l = 1), which gives the largest tachyonic mass, let us now consider the dependence on frequency 
in order to determine the full dispersion laws. With rj = cu/jkl one finds f38] 



a{ri) 



m 



{l + r]^ + e(-l + (6 + - (1 - 0^')) arctan 

T] + 1 + ie' 
1] — 1 + ie 



+ J^{r]'-1) l+er/'-(l + 0^1n 



Piv) 



2 2 

7] m 
2v^(l + er/T 



^1 + ^)(1 -^r]^) arctan 



1] + 1 + ie' 
7] — 1 + ie 



(17) 



The resulting poles in Ay and A^ are plotted in Fig. [6] by full and dashed lines, respectively, for the 
three cases ^ = (isotropic), ^ = —0.5 (prolate), and ^ = 5 (oblate). 

In the prolate case, ^ = —0.5, we see that there is magnetic screening down to and including the 
static limit, whereas the electrostatic screening is somewhat diminished (for fixed m). In the oblate 
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Figure 7: The growth rate 7 of the unstable modes in a Weibel instabihty as function of 
wave vector and for several values of anisotropy parameter ^. 



case, ^ = 5, the latter is instead enlarged but there is now a pole at zero frequency and positive k^, 
which is analytically connected with the dynamical magnetic screening poles at < by a line of 
poles with real k and negative uo"^. The corresponding imaginary values give the momentum- dependent 
growth rate 7(fc) of the unstable magnetic modes. These growth rates are shown in more detail in Fig. [7] 
for various values of ^. 

One can show that in the limit of large ^, the unstable modes are characterized by k^ax/m ~ i^^^-, 
but /c/m|^=^j^^^ ~ 1. Compared to the asymptotic gluon mass moo one has fcmax/'"^oo ~ and 

7maxMoo 1/^2- 



3 Discretized hard-loop effective theory 

In the Abelian case, the above dispersion laws give already a complete overview of the stable and 
unstable modes to leading order in the HL approximation. In particular, the unstable modes grow 
exponentially with rate 7(fc) until finally the limit of the HL approximation is reached and backreac- 
tion on the hard particle distribution /o has to be considered. However, with non- Abelian fields, the 
HL effective theory involves infinitely many vertex functions which become important when unstable 
modes have grown to amplitudes so large that the two terms in the covariant derivative 9^ — igA^ are 
comparable, i.e., when gA ~ m. This is still within the HL approximation, since hard particles have 
momentum m/g and so backreaction on them becomes important only when gA m/g ^ m. It is 
therefore an important question whether non-Abelian plasma instabilities are able to grow beyond the 
intrinsically non-Abelian regime gA ~ m, and this question is one that can be answered entirely within 
the HL approximation. 

In order to cope with the nonlocal structure of the HL effective theory, it is useful to introduce an 
auxiliary field formulation as developed in the HTL case in [HI |12] . Factorizing 

6r{x;p) = -^?iy;(t,x; v)9f^)/o(p) (19) 

one has 

[v ■ D{A)]W^{x; v) = F^^{A)v'' (20) 

and 

I^p(A)F- = ,-(-) = -9' I 7|^i^^^^^(-'-)- (21) 



9 



100 
10 

1 

0.1 
0.01 
0.001 
0.0001 

^^"°^0 2 4 6 8 10 12 14 

Figure 8: Lattice simulation of the non-Abelian Weibel instability with small initial seed 
fields that are constant in the transverse direction [15]. Plotted are the average energy 
densities £ in transverse/longitudinal chromomagnetic/electric fields and the total energy 
density transferred from hard particles to soft modes, £^(HL). 
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In the latter equation, only one linear combination of the four fields W^[x]\^ participates in the 
dynamical evolution. This extra field is a field on both, configuration space and velocity space, but in 
terms of this field together with the gauge potential, the dynamical equations are local and polynomial. 
These equations can therefore be solved by real-time lattice simulations upon discretizing both spacetime 
and velocity space. The induced current can then be written simply as 

3^{x) = ^Y^^^yV^ix) (22) 

by introducing M vectors v covering the unit sphere. In the HTL case, such a discretization has been 
studied first in |13] , whereas fl^ performed lattice simulations involving truncated series of spherical 
harmonics for that purpose. 

The first lattice study of non-Abelian plasma instabilities using discretized hard loops was performed 
in [IS] for modes with k || n. The restriction to modes which are constant with respect to spatial direc- 
tions perpendicular to n leads to the great simplification of a dimensional reduction to 1+1 spacetime 
dimensions that need to be discretized in addition to the 2-dimensional (compact) velocity space. The 
results obtained showed complicated nonlinear dynamics when Weibel instabilities enter the intrinsically 
non-Abelian regime, but except for a brief period of stagnation, it appeared that non-Abelian plasma 
instabilities continue to grow exponentially also in the late phase of the HL evolution, suggesting that 
they are essentially as effective for fast isotropization as their Abelian counterparts in conventional 
plasma physics. Closer inspection reveals that this continued growth is brought about by an effective 
Abelianization of the non-Abelian modes over extended (but finite) spatial domains. 

However, subsequent fully 3+1- dimensional lattice simulations [IHlll^ showed that this phenomenon 
is only true for the modes with exactly k || n. If there are also unstable modes with more general wave 
vectors, they eventually destroy the local Abelianization and lead to a saturation of the exponential 
growth, which goes over into a slow linear one (Fig. [2]). As Fig. ITU] shows, at the level of the individual 
modes this is accompanied by the build-up of an energy cascade where the energy fed into the soft 
unstable modes by the Weibel instability is distributed to higher-momentum stable modes through 
non-Abelian self-interactions [HI HQ] . 

More recent lattice simulations using larger lattices and larger (oblate) anisotropy Ref. |50| found 
a continued exponential growth of initially small perturbations in the case of very strong momentum 
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Figure 9: Lattice simulation of the non-Abelian Weibel instability as in Fig. [8] but now with 
generic, fully 3- dimensional initial seed fields and in linear scale, showing a saturation of 
exponential growth in the strong- field regime [47] . 
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Figure 10: Spectrum of field modes showing a saturation at low wave numbers and a cas- 
cading of energy to larger wave numbers with increasing simulation time (adapted from 
Ref. 08]). 



anisotropy, at least for small initial field configurations, similar to the behavior of the 1+1- dimensional 
case. With generic large initial gauge field amplitudes, again a linear growth was observed. 

4 Quark-gluon plasma instabilities in Bjorken expansion 

In the above treatment of anisotropic plasmas, a space-time independent background distribution was 
considered with a fixed momentum-space anisotropy. A somewhat more realistic case appropriate for 
the early stage of a quark-gluon plasma in formation after a collision of large nuclei is given by a 
distribution of hard massless particles which expands in one spatial dimension (the beam axis) in a 
boost-invariant manner. 

In the latter case it is convenient to switch to comoving spacetime coordinates, proper time r and 
rapidity 77, 

= (r = Vt^ - z^, x\ 7] = atanh(2/t)) (23) 
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with metric tensor Qa/s = diag(l, —1, —1, — t^). Momenta are parametrized by 

= |p_l| (cosh(?/ — r]), COS0, sin0, r~"^sinh(?/ — r])) (24) 



with y = atanh(p°/p^). A boost-invariant, transversely isotropic, free-streaming background distribu- 
tion is then given by 

/o(p, x) = /.o (v^i+P^) = /iso (sJpI + (P'^r/riso)^) (25) 



where p'^ is the boosted longitudinal momentum. Compared to f|T2l) we now have a proper-time depen- 
dent anisotropy parameter 

= (r/ri,o)' - 1, (26) 
such that there is an increasingly oblate momentum-space anisotropy for r > Tiso- This solves 

p^9^/o(p,x)=p"9„/o^ ^ ^ = (27) 

fixed 

because (p°'5a)p^(x) | fixed = 0. (Note that = — r^p''.) 

Clearly, we need to start at some finite tq > in order to avoid a singularity at r = 0. But a particle 
description of the initial stage of a heavy-ion collision can anyway make sense only after some finite time 
after the collision. In the so-called Color-Glass-Condensate (CGC) framework the earliest time when 
this is beginning to make sense is provided by the inverse saturation scale Q~^. In this framework, the 
initial hard-gluon density is given by [51j 

n(ro)=c (28) 

where c is the so-called gluon liberation factor. For definiteness, we shall assume Tq ~ in what 
follows, and we shall also assume that tq ^ Tiso, i.e., that the momentum distribution is always strongly 
anisotropic with positive ^. 

With increasing (proper) time, we have an increasing anisotropy with more and more modes becom- 
ing unstable, but at the same time the density decreases, which tends to diminish the growth rate. The 
competition between these effects will of course modify the evolution of non-Abelian plasma instabili- 
ties, but, at least in the free-streaming idealization, the HL framework can be generalized to study this 
case as well. 

Since p^di3[d'^p^fo{p±,prf)]\pt^=const. = (with index a upstairs!) we can again solve the Vlasov 
equation 

p-Z}5/,(p,x,t)^=,„„,,. =^/F|,a(-)/o(p,x,t), (29) 
by introducing auxiliary fields according to 

6nx;p) = -gW^iT,x\v;^,y)d^p)fo{p±,Pv) (30) 

with satisfying 

V ■ D Wo^{t, x\ ri- 0, y%,y = v^F^p, (31) 

where we now define 

" - ^" ( \.( \ A ■ - sinh(?/-r7)x 
V = - r = (cosh (y — r]), cos 0, sm0, . (61) 

|P±I ^ T ^ 

The induced current in the non-Abelian Maxwell equations, which in comoving coordinates now 
read 

-D^{tF^^)=3^ (33) 
r 

12 



20 40 60 80 100 

X/X„ 



Figure 11: Results [54] from a 1D+3V real-time lattice simulation of non-Abelian plasma in- 
stabilities in Bjorken expansion, seeded by small initial rapidity fluctuations with a spectrum 
modelled after Ref. [55]. The plot shows the proper-time dependence of the total chromo- 
field energy density and its individual components £ = Sbt + ^Et + ^Bl + ^El = £t + £l as 
well as the chromo- field energy gain rate R defined hj R = dS/dr + 2St/t. 



is given by 

-D^{TF'^P)=f{T,x\v) = r dyv^W{T,x\v-<P,y) (34) 

r Z Jo zir J-oo 

with 

W = v'Wi - 4- sinh(?/ -'n)W^. (35) 

''"iso 

In the expanding case, already the effectively Abelian, weak-field behavior is highly nontrivial, 
leading to complicated integro-differential equations [52], from which one can extract that Weibel in- 
stabilities behave asymptotically like 

I*(r,z/) ~r 2F3(^,^;2,2-zz/,2 + iz/;-/ir) , s = Vl + 4z/2 (36) 

where v is the wave number corresponding to the rapidity variable ?7, and 

= ^Tisom^. (37) 

o 

For both v and r large, this leads to 

I^(r) ~ r^/^exp(2yA^), (38) 

which can be understood by the large-^ limit of 7 quoted at the end of Sect. 12.31 together with the 
fact that the thermal mass scale drops like the square root of hard particle density, which behaves 
as n ~ 1/r. Such a behavior was found previously in numerical CGC simulations that include small 
rapidity fluctuations [53] . 

In order to also investigate specific non-Abelian dynamics in the HL framework, it is again necessary 
to discretize space-time, now parametrized by proper-time and rapidity coordinates, and the non- 
compact velocity/momentum-rapidity space parametrized by and y. In fact, because the integrand in 
f lMl) is exponentially suppressed at large |?/ — ?7|, only a finite rapidity interval is required in numerical 
simulations. 

Fig. E] shows the results of a real-time lattice simulation [54j of the expanding HL equations (I3T!) 
and (134|1 with small initial seed fields that are constant in transverse space and supplied with an initial 
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spectrum inspired by the CGC scenario [55] and a gluon-liberation factor [56] c = 2 In 2 (which is also 
close to the recent numerical results of Ref. [57j). Given the findings in stationary anisotropic plasmas, 
the growth of non-Abelian plasma instabilities in the 1+1-dimensional setting probably gives an upper 
bound on more generic cases. Considering that for — 1 GeV for RHIC and ~ 3 GeV for the LHC 
and that thus the maximal time in Fig. [TT] corresponds to some 20 fm/c for RHIC and 7 fm/c for the 
LHC, one finds an uncomfortably long delay for the onset of plasma instabilities at least for RHIC. 
Further studies of more generic initial conditions (including then also strong initial fields) are clearly 
called for to investigate this situation in more generality. 

5 Conclusion 

The extension of the HTL effective theory to anisotropic plasmas leads to a (leading-order) theory of 
non-Abelian plasma instabilities, which in the weak-field limit can be analysed by the study of dispersion 
laws, but which in the strong-field case requires numerical simulations on real-time lattices. The latter 
have shown that contrary to initial expectations non-Abelian plasma instabilities tend to saturate in 
the nonlinear regime where fields are nonperturbatively large, but not yet large enough for immediately 
modifying the anisotropic background distribution of hard particles. However, a complete isotropization 
scenario still needs to be worked out. Moreover, when boost-invariant expansion is taken into account, 
there is an uncomfortably long delay for the onset of growth, which seems too long for the environment 
provided by RHIC, though not necessarily for the LHC, which may be where wQGP physics eventually 
comes into its own. 

It should be kept in mind that the above studies are based on a truly weakly coupled situation with 
g ^ 1 and it is perhaps not so surprising that at least RHIC physics lies outside of the results obtained 
by extrapolation to g ^ 1, since this is also the case when perturbative thermodynamic bulk results 
are extrapolated to below 2-3 Tc. Indeed, numerical simulations which do not separate collective HL 
physics from other interactions have concluded that fast thermalization through plasma instabilities 
may be indeed possible, see [55] . 

Still, it remains a theoretical challenge to understand these issues in the limit (7^1 where a system- 
atic perturbative analysis should be possible, and, as said, perhaps this turns out to find applications 
eventually at the higher scales to be probed by the LHC. 

On the theoretical side, a systematic treatment of the HL effective theory also offers new concep- 
tual problems, as shown by its application to heavy-ion observables such as energy loss [59] and the 
momentum broadening of jets [60l [61], [62]. There is clearly more work to be done, both numerically 
and analytically. 
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